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In the present work we investigate one possible variation on the usual electrically neutral pulsars: 
the inclusion of electric charge. We study the effect of electric charge in pulsars assuming that the 
charge distribution is proportional to the energy density. All calculations were performed for zero 
temperature and fixed entropy equations of state. 
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I. INTRODUCTION 

Pulsars are believed to be the remnants of supernova 
explosions. They have masses 1-2Mq, radii ~ 10 km, 
and a temperature of the order of 10 11 K at birth, cooling 
within a few days to about 10 10 K by emitting neutrinos. 
Pulsars are normally known as neutron stars. Qualita- 
tively, a neutron star is analogous to a white dwarf star, 
with the pressure due to degenerate neutrons rather than 
degenerate electrons. The assumption that the neutrons 
in a neutron star can be treated as an ideal gas is not well 
justified: the effect of the strong force needs to be taken 
into account by replacing the equation of state (EoS) for 
an ideal gas by a more realistic EoS. The composition 
of pulsars remains a source of speculation, with some of 
the possibilities being the presence of hyperons [ll [2, 0] , 
a mixed phase of hyperons and quarks 0, [E lEMHii a 
phase of deconfmed quarks or pion and kaon condensates 
[9j. Another possibility would be that pulsars are, in fact, 
quark stars [10] . In conventional models, hadrons are as- 
sumed to be the true ground state of the strong interac- 
tion. However, it has been argued (HI [I2> EE 13 EE] that 
strange matter composed of deconfmed u, d and s quarks 
is the true ground state of all matter. In the stellar mod- 
eling, the structure of the star depends on the assumed 
EoS, which is different in each of the above mentioned 
cases. An important distinction between quark stars and 
conventional neutron stars is that the quark stars are self- 
bound by the strong interaction, whereas neutron stars 
are bound by gravity. 

Once an adequate EoS is chosen, it is used as in- 
put to the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions [l6j], which are derived from Einstein's equations 
in the Schwarzschild metric for a static, spherical star. 
Some of the stellar properties, as the radius, gravitational 
and baryonic masses, central energy densities, etc are ob- 
tained. These results are then tested against some of the 
constraints provided by astronomers and astrophysicists 
[i~7l [Lsj and some of the EoS are shown to be inappropri- 
ate for describing pulsars [E IE HI • 

Notice also that the temperature in the interior of the 
star is not constant |8|,l21(, but the entropy per baryon is. 
This is the reason for choosing fixed entropies to take the 
temperature effects into account. The maximum entropy 



per baryon (S) reached in the core of a new born star is 
about 2 (in units of Boltzmann's constant) [IE]. We then 
use EoS obtained with S = (T = 0), 1 and 2. 

We study the effects of the electric charge in compact 
stars. This study was first performed in stars composed 
of hot ionized gas [12] and then reconsidered for a cold 
star (T = 0) described by a polytropic EoS [H. The 
electric charge distribution is assumed to be proportional 
to the mass density. The TOV equations again have to 
be modified to take the electric field into account. In the 
present work once more all possible classes of pulsars are 
examined in the presence of the electric field for S — 0, 1 
and 2. 

This paper is organized as follows: in Sec. II the for- 
malisms of the electrically charged stars are revisited and 
the results are presented. In Sec. Ill the results are dis- 
cussed and the main conclusions are drawn. 



II. FORMALISM AND RESULTS 

As the first step we need to know the EoS of the system, 

e=e(p), n = n(p), 

where p is the pressure, e is the energy density, and n is 
the number density of baryons. Once an adequate EoS is 
obtained, it can be used to provide the stellar properties. 

A. Electrically charged compact stars 

In this section we include modifications in the TOV 
equations to describe electrically charged pulsars with 
null angular velocity. The geometry that describes a 
static spherical star is given by equation (1). In order 
that the Maxwell equations are incorporated into the 
stress tensor it becomes: 



Tt ' = { P + e)u»u„-p6: 

w 

where again p is the pressure, e is the energy density, and 
is the 4- velocity vector. 
The electromagnetic field obeys the relation 



-gF^} 



(2) 



where is the four current density. Next we consider 
static stars only. Hence the electromagnetic field is only 



2 



due to the electric charge, which means that F 01 = — F 10 , 
and the other terms are absent. From the four-potential 
the surviving potential is Aq = <\>. Thus the electric 
field is given by 



E( r ) = 1 [\irr 2 fe^ +x V 2 dr, 
r Jo 



(3) 



where f'e v l 2 = p c h is the charge density. The electric 
field can be written as 



dE(r) 2E A/2 

— ; — = 1- 47rp c/l e ' 

dr r 

and total charge of the system as 



Q 



4nr 2 p ch e x/2 dr, 



(4) 



(5) 



where R is the radius of the star. 
In the star frame the mass is 

dM tot 2 / E(r) 2 

__( r ) = 47rr [e+ — 

To an observer at infinity, the mass is 
E(r) 2 



(6) 



Moo = / 47rr 
'0 



8tt 



dr = Mtot(R) 



Q{R? 

2R 



(7) 

By using the conservation law of the stress tensor 
(T^ M = 0) we obtain the hydrostatic equation 



dp 
dr 



Airr ( p 



E(rY 
8tt 



P) 



r 2 (I 



2M to 



+ p ch E{r)e x / 2 . 



(8) 



The first term on the right-hand side comes from the 
gravitational force and the second term comes from the 
Coulomb force. By using the metric and the relation 



1 



K - 5^ 



we obtain the following differential equation 
E{r) 2 



dr 



8nre 



1 



(9) 



(10) 



which is used to determine the metric e x . 

So, we have a set of differential equations to be solved 
formed by equations(19), (21), (23) and (25). The bound- 
ary conditions at r — are E(r) — 0, e x = 1, n = p c and 
at r = R, p = 0. We assume that the charge goes with 
the energy density e as prescribed in [20] : 



Pch = / x 0.86924 x 10 3 e. 



(11) 



This choice of charge distribution is a reasonable as- 
sumption in the sense that a large mass can hold a 



large amount of charge. In table Q] results f° r electrically 
charged neutron stars are presented. We have calculated 
the results for 48 different configuration models of com- 
pact stars. Once again the EoS for hadronic and hybrid 
stars were taken from Q, the EoS for quarkionic stars 
were taken from [l(| . In table Q] the electric charge Q is 
given in Coulomb and / varies from zero (no charge) to 
a small value (0.0006). The related mass-radius plots for 
hadronic, hybrid and quarkionic stars are given respec- 
tively in figs. [HU and EE] 



III. RESULTS AND CONCLUSIONS 

Let's now go back to our results in order to compare 
them with what is found in the literature and draw the 
conclusions. 

In [23J the models named 0,P,Q and R given in table 
[J can again be compared with our result for the hadronic 
star at T = and they are indeed very similar. 

We next look at the results displayed in table Q] and 
corresponding figures. The general trend is the same ob- 
served in a simple polytropic EoS for T = [20| . i.e., 
the electric charge, the maximum mass and the mass ob- 
served at infinity increase with /, as it should be. Al- 
though the EoS used in the present work are very dif- 
ferent from the one used in [2(| the values of the radii 
obtained and the electric charge for a fixed / value are 
compatible. Figs. IH2l3l andl4lalso show the same behav- 
ior as fig. 2 of [20J, i.e., as / increases, the maximum mass 
and radius of a family of stars increase. The EoS used 
in the present work are for bare pulsars, i.e. the outmost 
layer is not included. From table [J one can see that the 
effect of entropy on a charged star remains the same as 
in a neutral star: the maximum masses and the radii de- 
crease with the increase of the entropy for hadronic and 
quarkionic stars within the NJL model. For hybrid and 
MIT stars the behavior is not so well defined. 

In the present work we have investigated one possible 
variation on the usual electrically neutral pulsars: the 
inclusion of electric charge. We have observed that the 
behaviors shown in previous works with much simpler 
EoS were also observed here. The influence of the tem- 
perature was also investigated. We are now in a position 
to calculate the energy released from the conversion of 
a metastable star (hadronic or hybrid) to a stable star 
(hybrid or quarkionic) under the influence of the electric 
charge. A more detailed and complete investigation, with 
a smaller bag parameter in the MIT bag model is under 
way. 
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FIG. 3: Solutions for electrically charged quarkionic stars obtained with the MIT bag model for different values of /. 
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TABLE I: Electrically compact stars with different charge fraction /. 



Type 


Entropy 


/ 


M max 


Moo 


R 


e c 


Q 








(Mo) 


(M ) 


(km) 


(g/cm 3 ) 


(C) 


Hadronic 








2.04 


2.04 


11.72 


1.98 x 10 15 





Hadronic 





0.0002 


2.08 


2.10 


11.84 


1.94 x 10 15 


7.97 x 10 19 


Hadronic 





0.0004 


2.22 


2.28 


12.18 


1.84 x 10 15 


1.71 x 10 20 


Hadronic 





0.0006 


2.50 


2.66 


12.73 


1.75 x 10 15 


2.91 x 10 20 


Hadronic 


1 





f.96 


1.96 


11.02 


2.23 x 10 15 





Hadronic 


1 


0.0002 


2.00 


2.02 


11.15 


2.13 x 10 15 


7.69 x 10 19 


Hadronic 


1 


0.0004 


2.13 


2.19 


11.44 


2.04 x 10 15 


1.64 x 10 20 


Hadronic 


1 


0.0006 


2.39 


2.55 


12.01 


1.85 x 10 15 


2.78 x 10 20 


Hadronic 


2 





1.93 


1.93 


10.91 


2.24 x 10 15 





Hadronic 


2 


0.0002 


1.97 


1.98 


11.01 


2.19 x 10 15 


7.55 x 10 19 


Hadronic 


2 


0.0004 


2.09 


2.15 


11.26 


2.15 x 10 15 


1.61 x 10 2() 


Hadronic 


2 


0.0006 


2.34 


2.50 


11.86 


1.90 x 10 15 


2.72 x 10 20 


hybrid 








1.64 


1.64 


12.33 


1.57 x 10 15 





hybrid 





0.0002 


1.68 


1.69 


12.43 


1.57 x 10 15 


5.98 x 10 19 


hybrid 





0.0004 


1.82 


1.86 


12.82 


1.48 x 10 15 


1.31 x 10 20 


hybrid 





0.0006 


2.13 


2.23 


13.52 


1.39 x 10 15 


2.31 x 10 20 


hybrid 


1 





1.50 


1.50 


11.32 


1.75 x 10 15 





hybrid 


f 


0.0002 


1.54 


1.55 


11.43 


1.71 x 10 15 


5.44 x 10 19 


hybrid 


f 


0.0004 


1.67 


1.70 


11.74 


1.66 x 10 15 


1.19 x 10 20 


hybrid 


f 


0.0006 


1.94 


2.03 


12.34 


1.58 x 10 14 


2.10 x 10 20 


hybrid 


2 





1.50 


1.50 


11.76 


1.58 x 10 15 





hybrid 


2 


0.0002 


1.54 


1.55 


11.86 


1.58 x 10 15 


5.41 x 10 19 


hybrid 


2 


0.0004 


1.68 


1.71 


12.21 


1.53 x 10 15 


1.18 x 10 20 


hybrid 


2 


0.0006 


1.95 


2.04 


12.86 


1.44 x 10 14 


2.10 x 10 20 


Quarkonic(MIT) 








1.22 


1.22 


6.77 


5.14 x 10 15 





Quarkonic(MIT) 





0.0002 


1.25 


1.26 


6.81 


5.13 x 10 1B 


4.75 x 10 19 


Quarkonic(MlT) 





0.0004 


1.33 


1.37 


6.97 


4.95 x 10 15 


1.02 x 10 20 


Quarkonic(MIT) 





0.0006 


1.50 


1.60 


7.28 


4.56 x 10 15 


1.74 x 10 20 


Quarkonic(MIT) 


1 





1.22 


1.22 


6.76 


5.17 x 10 15 





Quarkonic(MIT) 


f 


0.0002 


1.25 


1.26 


6.82 


5.07 x 10 15 


4.75 x 10 19 


Quarkonic(MlT) 


1 


0.0004 


1.33 


1.37 


6.98 


4.88 x 10 15 


1.02 x 10 20 


Quarkonic(MlT) 


1 


0.0006 


1.50 


1.60 


7.29 


4.50 x 10 1B 


1.74 x 10 20 


Quarkonic(MIT) 


2 





1.23 


1.23 


6.79 


5.08 x 10 15 





Quarkonic(MlT) 


2 


0.0002 


1.25 


1.26 


6.83 


5.09 x 10 15 


4.77 x 10 19 


Quarkonic(MlT) 


2 


0.0004 


1.34 


1.37 


6.98 


4.94 x 10 15 


1.03 x 10 20 


Quarkonic(MfT) 


2 


0.0006 


1.50 


1.60 


7.31 


4.49 x 10 15 


1.74 x 10 20 


Quarkonic(NJL) 








1.20 


1.20 


7.87 


3.45 x 10 15 





Quarkonic(NJL) 





0.0002 


1.23 


1.23 


7.93 


3.45 x 10 15 


4.42 x 10 19 


Quarkonic(NJL) 





0.0004 


1.32 


1.34 


8.14 


3.30 x 10 15 


9.54 x 10 19 


Quarkonic(NJL) 





0.0006 


1.49 


1.57 


8.48 


3.20 x 10 15 


1.64 x 10 20 


Quarkonic(NJL) 


1 





1.17 


1.17 


7.71 


3.68 x 10 15 





Quarkonic(NJL) 


1 


0.0002 


1.19 


1.20 


7.79 


3.56 x 10 15 


4.32 x 10 19 


Quarkonic(NJL) 


1 


0.0004 


1.28 


1.31 


8.00 


3.50 x 10 1B 


9.31 x 10 19 


Quarkonic(NJL) 


f 


0.0006 


1.45 


1.53 


8.31 


3.38 x 10 15 


1.60 x 10 20 


Quarkonic(NJL) 


2 





1.10 


1.10 


7.18 


4.58 x 10 15 





Quarkonic(NJL) 


2 


0.0002 


1.13 


1.13 


7.29 


4.34 x 10 15 


4.09 x 10 19 


Quarkonic(NJL) 


2 


0.0004 


1.20 


1.23 


7.42 


4.40 x 10 15 


8.23 x 10 19 


Quarkonic(NJL) 


2 


0.0006 


1.36 


1.43 


7.80 


3.98 x 10 15 


1.50 x 10 20 



